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Abstract: We present fully self-consistent calculations in the 
fluctuation exchange approximation for the half-filled Hubbard 
model in 2D. A non-fermi liquid state evolves with decreasing 
temperature in this self-consistent model of coupled spin fluc- 
tuations and quasiparticles. The mean field phase transition 
to long-range antiferromagnetic order is suppressed and we find 
no evidence of a phase transition to long-range magnetic order. 
We show that the real part of the self-energy at zero energy 
shows a positive slope and the imaginary part of the self-energy 
shows a local minimum. The scale of this structure is set by 
the zero temperature gap in mean field theory. The growth of 
spin fluctuations is reflected in the evolution of sharply peaked 
structure in the spin fluctuation propagator around zero energy 
and Q = (vr, vr). We present calculations for the Hubbard model 
in ID in a two-particle self-consistent parquet approximation. 
A second moment sum rule suggests that vertex corrections to 
the self-energy are responsible for the greater accuracy of the 
parquet approximation as compared to other self-consistent per- 
turbation theories. 

§ 1 Introduction and Physical Picture 

The role of spin fluctuations in the physics of a many-body Fermi system seems a 
particularly appropriate topic for a festschrift in honor of David Pines. In liquid ^He, 
ferromagnetic spin fluctuations play a central role in the normal and superfluid states. 
The exchange of spin fluctuations contributes a non-analytic correction to the quasipar- 
ticle energy that is reflected in finite-temperature corrections to the transport properties 
and specific heat predicted by Landau Fermi liquid theory (Baym and Pethick, 1991). 
Similar corrections were also found for metallic systems near a ferromagnetic instability 
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(Berk and Schrieffer 1966, Doniach and Englesberg 1966). In both heavy fermion materi- 
als and the high-tempcraturc superconductors, the close proximity of antiferromagnetic 
instabilities and the observation of unusual and possibly unconventional superconduc- 
tivity has aroused interest in the role of antiferromagnetic spin fluctuations, both as an 
explanation of the unusual normal state properties and as a pairing mechanism. 

The normal-state properties of the high-temperature superconductors have inspired 
speculations that interplay among correlations, low dimensionality and the tendency to 
magnetic order could give rise to low-energy excitations that are not adequately de- 
scribed by the Landau theory of a Fermi liquid. Varying degrees of departure from 
Landau theory have been proposed (for recent reviews, see Kampf (1994) and Dagatto 
(1994)), ranging from the 'soft singularities' of phenomenological marginal Fermi liquids 
and nearly antiferromagnetic Fermi liquids to more dramatic conjectures of separation 
of charge and spin excitations such as Luttinger liquids (Anderson, 1990). All the pro- 
posed departures from a Fermi liquid description of low energy excitations ultimately 
rest on subtle features of interactions between renormalized quasiparticles. Thus a self- 
consistent approach seems necessary for a complete and satisfactory theory. The obser- 
vation of an antiferromagnetically ordered state in the parent compounds of the high 
temperature superconductors suggests the possibility that the renormalization of quasi- 
particle excitations by the exchange of spin fluctuations may account for the unusual 
properties of these materials. Here we describe results for a self-consistent model aimed 
at elucidating the physics of coupled spin fluctuations and quasiparticle excitations. 

The 2D copper-oxide layer structure of the cuprate high temperature superconduc- 
tors has inspired great interest in the single-band Hubbard Hamiltonian in 2D as a model 
for the high temperature superconductors (Anderson, 1987). Numerical evidence from 
quantum Monte Carlo calculations strongly indicates that the ground state of the half- 
filled Hubbard model in 2D possess long-range antiferromagnetic order (Hirsch and Tang, 
1989; White, Scalapino, Sugar, Loh, Gubernatis, and Scalettar, 1989); at finite temper- 
ature, the possibility of long-range magnetic order is precluded by the Mermin- Wagner 
theorem. In this paper we consider the 2D Hubbard model at half-filling, where spin 
fluctuations are expected to play a pivotal role as the system of coupled quasiparticles 
and (incoherent) spin fluctuations self-consistently evolves with decreasing temperature 
toward the ordered state at zero temperature. For small U. mean field theory predicts 
a phase transition to an antiferromagnetically ordered state at a non-zero temperature. 

In a previous work (Deisz, Hess and Serene, 1995) we showed that the self-consistent 
inclusion of fiuctuations at the level of the fiuctuation exchange approximation (FEA) 
strongly suppresses the mean field phase transition; we have not observed a magnetic 
phase transition down to the lowest temperatures for which we have calculations. At 
sufficiently low temperature, the coupling of growing spin fluctuations to quasiparticle 
excitations leads to the evolution of a non-Fermi liquid state, clearly signaled by anoma- 
lous features in the self-energy for £ ~ 0: a positive slope in the real part and a local 
minimum in the imaginary part. Spectral weight is rapidly pushed away from the Fermi 
surface with decreasing temperature leading to the formation of a weak pseudogap in the 
single-particle density of states. We presented a simple analytical model which showed 
that these features, along with developing shadow structure in single-particle spectral 
functions, are a consequence of evolving spin fluctuations and their self-consistent effect 
on the structure of evolving quasiparticles. 

This paper extends our previous results providing further details of the evolution of 
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the non-Fermi liquid state with decreasing temperature. We observe that the structure 
that develops in the self-energy at low energy occurs on the scale of the zero temperature 
gap in mean field theory. The formation of the non-Fermi liquid state is accompanied 
by the evolution of a sharply peaked structure in the T-matrix; the temperature and 
energy dependence of this structure suggests the existence of a strongly temperature 
dependent characteristic energy scale not unlike that proposed in the phenomenologically 
motivated NAFL (Pines, 1990; Millis, Monein and Pines, 1990). To make this discussion 
self contained, we briefly summarize the FEA in the next section. In §3 we present a 
short summary of our numerical methods; it will be convenient there to summarize some 
formally exact results that suggest that FEA is a better approximation than including 
only the sum of the particle- hole bubble diagrams (the shielded potential approximation). 
In §4 we describe the non-Fermi liquid solution that evolves at low temperature. Finally 
in §5, we present results for the ID Hubbard model at half filling in a fully self-consistent 
parquet approximation. This approximation includes interactions between fluctuations 
in the calculation of the self-energy and is self-consistent at the level of the two-particle 
propagator. 

§^ The Fluctuation Exchange Approximation for the Hubbard Model 

In the propagator renormalized perturbation theory of Luttinger and Ward (Lut- 
tinger and Ward, 1960) the grand thermodynamic potential Q is constructed as an 
independent functional of the fully renormalized single particle propagator G and the 
self-energy E, and is stationary with respect to independent variations of G and E. 
Dyson's equation follows from 

^^^0. ,1) 

and the representation of the self-energy as a sum of 'skeleton' graphs evaluated with 
the fully renormalized propagator and bare two-particle vertex is obtained from 

^ ^ 0. ,2) 

The self-consistent solution of Eq. (|l]) and Eq. (|) yields the exact Green's function. 
Explicit summation of the infinite set of 'skeleton' graphs for the self-energy required 
for the exact solution is generally not possible and approximations are generated by 
selecting a class or classes of vacuum graphs for the grand potential. The FEA is 
one such approximation; it includes the exchange of density fluctuations, spin-density 
fluctuations and singlet-pair fluctuations. 

It is convenient to discuss the FEA by alternating between position and (imaginary) 
time, and momentum and (imaginary) frequency representations. The FEA for the 
electron self-energy of the paramagnetic state is given by 

S(r, r) = f/2 [xphi^, r) + T,,(r, r) + T,/(r, r)] ^(r, r) + U%p{r, r)G'(-r, -r) (3) 

where Xph = — G'(r, r) (?(— r, — r) and Xpp = G{r,T) G{r,T) are the particle-hole and 
particle-particle susceptibility bubbles. The density, spin, and particle fluctuation T- 
matrices are most simply written in momentum-frequency space, 

\2 



Tsf{q,uJm) = - r 4 

2 l-UXph{(l,i^m) 
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Tpp{q_,uj^) = -- — — (5) 

2 1 + UXph{(l,i^m) 

rr ( \ UXpp{fl,^mf 

= l + t/x>(q.^„) - 

The Green's function is obtained from Dyson's equation 

^(k, = G„(k, - E(k, £„) (7) 

where G'o(k, £„) is the non-interacting Green's function. The FEA self-energy (and 
therefore G) are found from the self-consistent solution of Eqs. (§| - 0). 

§5 Numerical Methods and Analytic Considerations 

We briefly describe our methods here; a more detailed account can be found else- 
where (Deisz, Hess and Serene, 1994). The self-consistent self-energy is obtained by 
the iterative solution of Eqs. (0-0). Efficient fast Fourier transforms are used to switch 
between (r, r)-space and (k, e„)-space. Since we use periodic boundary conditions on 
a periodic lattice, the use of discrete Fourier transform (DFT) methods introduces no 
fundamental difficulties in the transformation between k and r variables. On the other 
hand, a discrete sampling of the r variable generally introduces an artificial periodicity 
(with a period of the cutoff frequency) in the space of (discrete) Matsubara frequency. 
Nonetheless, all currently used high frequency cutoff schemes can be simulated using 
DFT's, with the possible cost of doubling the number of frequencies used (Serene and 
Hess, 1992). Building on the observation that comparatively little information is con- 
tained in the high frequency behavior of G, the T-matrices, the susceptibility bubbles, 
and S, we have devised a substantially more accurate method that preserves the formal 
internal self-consistency of the theory, improves the stability of the algorithm, and en- 
ables thermodynamic quantities to be calculated rehably. We decompose all quantities 
into the sum of a numerical term and an analytic term. The analytic terms are asymp- 
totically exact at high frequency; they allow highly accurate calculations of contributions 
to low-frequency properties from convolution sums over high-frequency tails. These con- 
tributions, which propagate through any iterative solution of the nonlinear equations of 
the FEA, are lost in other implementations that introduce a high frequency cutoff. We 
find that accurate treatment of these contributions is absolutely essential for meaningful 
results from the FEA in the presence of large fiuctuations. We perform calculations on 
128 X 128 lattices of k-points with typically 512 Matsubara frequencies. A fine mesh is 
required to represent accurately important fine structure that evolves, especially in the 
spin-fiuct nation propagator, with decreasing temperature. The resulting algorithm has 
yielded converged solutions at half-filling for a large range of U with continuous control 
of the temperature T down to T ~ 0.008t the lowest temperatures we have explored. 
Numerical analytic continuation of quantities calculated on the imaginary frequency 
axis to the real frequency axis is accomplished using Fade approximants (Vidberg and 
Serene, 1977). 

The slowly decaying high frequency tails of the single-particle Green's function are 
a consequence of the discontinuity of G and its derivatives at r = 0. For example, the 
discontinuity G{0'^) — G{0~) = —1 ( which follows from the canonical anticommutation 
relations for Fermions) results in the l/i^n behavior of G at high frequency. More 



4 





Figure 1: The evolution of a non- Fermi liquid self-energy with decreasing 

temperature. Shown arc S(kp' = (vr/S, Tvr/S), e) for T = 0.06 (long dash), 
0.07 (dashed), 0.09 (sohd) and 0.12 (dot-dash). Note that above e ~ 1, Im 
S(k, e) displays the linear behavior consistent with our earlier observations 
(Serene and Hess, 1992). 

generally, discontinuities in G{t) and its derivatives can be related to weighted integrals 
of the spectral function, 

- q^^|„„- = (-!)-'> X>M(M)<<.. (8) 

The left hand side of this equation can be evaluated directly using the equations of motion 
for G, yielding exact sum rules for the single particle spectral function. From Dyson's 
equation for the exact G and the 'skeleton' diagram expansion for S one can show (Dcisz, 
Serene and Hess, 1995) that for the Hubbard Hamiltonian: (1) the contribution from the 
self-energy to the n = 1 sum rule comes entirely from the Hartree-Fock diagram with 
the exact G, and (2) the contribution from the self-energy to the n — 2 sum rule comes 
entirely from the second-order 'skeleton' diagram evaluated with the exact G. When 
evaluated with the FEA Green's function, the first moment sum rule is satisfied but 
the second generally is not. It is important to note, however, that the violation of the 
sum rule is significantly greater for approximations that keep only particle-hole bubble 
contributions. From this point of view, the FEA is a better approximation than taking 
particle-hole bubble or particle-particle ladder 'skeleton' graphs alone. 

^4 Neither Fermi-Liquid Nor Antiferromagnet 
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Figure 2: Im T,(kF,e) for U — 2.7 and T = 0.10 showing the smooth trans- 
formation from clear non-fermi hquid behavior at k = (vr, 0) on the Fermi 
surface (lowest curve at £ = 0) to a Fermi-liquid like behavior at the point 
k = (7r/2,7r/2) on the Fermi surface along the (1, 1) direction (highest curve 
at £ = 0). 

We have calculated self-consistent self-energies in the fluctuation exchange approx- 
imation for the 2D Hubbard model at half filling. In Fig. 1, we show the self-energy 
at low energy for a momentum on the Fermi surface near the X point and for various 
temperatures. For this modest U = 1.57 (all energies are understood to be measured 
with respect to the hopping matrix element t) we observe that at high temperature 
(T ~ 0.12, approximately the mean field transition temperature) the energy dependence 
of E at low energy is roughly similar to that of Fermi hquid theory: the slope of Re 
E corresponds to a quasiparticle pole weight near unity, and Im E shows a roughly 
parabolic energy dependence. It is important to note, however, that Im S(k^,£ = 0) is 
significantly larger in magnitude compared with that calculated for an approximation 
containing only the second order diagram (Deisz et ai, 1995). This refiects enhanced 
quasiparticle-quasiparticle scattering by AFM spin fluctuations. As the temperature is 
lowered, E develops structure at £ = that is inconsistent with a Fermi liquid; the 
slope of Re S decreases rapidly and eventually changes sign, and Im S develops a lo- 
cal minimum. These anomalies occur at roughly half the mean field theory transition 
temperature to the antiferromagnetically ordered state in mean field theory. Moreover, 
the energy scale suggested by the local maxima in Im E at £ ~ 0.28 is very nearly the 
size of the zero temperature gap A in mean field theory obtained by the k-sum over the 
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magnetic Brillion zone, 



u- 



N 



E 



(9) 



Wc have found this near agreement to hold for several values of U less than half of the 
bandwidth. In a separate work, (Deisz,et. ai, 1995) we have shown from calculated 
single-particle spectral functions that the evolution of the anomaly is accompanied by a 
dramatic suppression of spectral weight at the Fermi surface. The loss in spectral weight 
at the Fermi surface leads to a formation of a weak pseudogap in the density of states 
at e = 0. As shown in Fig. 2, the evolution of the anomaly is anisotropic; it starts at 
the X point and spreads across the Fermi surface with decreasing temperature. 
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Figure 3: The monotonic evolution with decreasing temperature of sharply 
peaked structure in the spin-fluctuation T-matrix Tgf on the real energy 
axis at Q = {11,71). Shown are Re Tsf{Q,uj) and Im Ts/(Q,ci;) for T — 
0.20, 0.18, 0.10, 0.08, 0.07, 0.065, 0.06 

The appearance of the non-Fermi liquid state is presaged by the evolution and rapid 
growth with decreasing temperature of structure in the spin-fluctuation T-matrix that is 
sharply peaked in momentum at Q = (tt, tt) and in Matsubara frequency about Um = 0. 
The analytic continuation of Tsf{Q,ujm) to the real frequency axis is shown in Fig. 3. 
The rapid increase in Tsf{Q,uj = 0) is shown in detail together with the temperature 
dependence of the width of the peak in frequency A^T^j and momentum AqT^j (for high 
temperature, the latter quantity is anisotropic and the width is for the kx direction). 
Also plotted in the lower panel of Fig. 4 is l/Tsf{Q,u! = 0). Below T ~ 0.07, A^^Tgf 
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and l/Tsf{Q,u! = 0) track each other closely suggesting the existance of a characteristic 
temperature T* below which the behavior of Tsf{Q, u) can be characterized by a simple 
diffusive pole at a (temperature dependent) characteristic energy lj*. It is interesting to 
note that uj* decreases rapidly with decreasing temperature below T* and the admittedly 
small amount of data available below T* can be approximated well by the functional 
form Tsf{Q,0) = B exp{A/T). We note that the structure of Tsj(q ~ Q,^^) below 
T* is similar to that of the spin-spin correlation function proposed by Millis, Monein 
and Pines, 1989. In contrast, Tsf{q,uj) presented here is the result of a fully self- 
consistent calculation and T's/(q, u) cannot be observed directly in an experiment; a fully 
self-consistent calculation of the dynamical spin-spin response function is necessary to 
determine how u* is expressed in measurable magnetic properties. 

We emphasize that the rapid rise in Tsf{Q, 0) with decreasing temperature need not 
be reflected in any physically measurable response function and does not necessarily 
signal a phase transition. The fully self-consistent staggered spin susceptibility as a 
function of temperature is presented in Fig. 5. The staggered spin susceptibility was 
calculated from a trace of the self-consistent Green's functions in the presence of an ap- 
plied staggered field. The absence in Xafm of any reflection of the rapid rise in T^f^Q, 0), 
dramatically contrary to expectations from an RPA calculation of the response, shows 
the importance of vertex corrections in the FEA response functions near an instability. 

§5 Beyond FEA: Calculations in Parquet Approximations 

De Dominicis (De Dominicis, 1963) reformulated the propagator renormalized pertur- 
bation theory of Luttinger and Ward to eliminate not only the single-particle potentials 
in favor of single-particle Green's function but also to eliminate the two-body poten- 
tial in favor of the two-particle propagator G2 (or equivalently the two-particle vertex 
function). In this theory, the entropy S plays a role analogous to that of the grand 
potential in the theory of Luttinger and Ward; the entropy may be viewed as a func- 
tional of G, S, G2 and the two-particle vertex function F. The stationarity conditions 
for 5* with respect to variations of G, S, G2 and F constitute a closed set of equations 
that must be solved self-consistently. Here we present the equations for the two-particle 
self-consistent perturbation theory schematically and refer the interested reader to the 
original works for the full formulation. These are the: Dyson's equation for the single 
particle propagator in Eq. (|^); a relation between the self-energy, the bare interaction 
U, and the renormalized two-particle vertex function 

E = UGGGT, (10) 

which includes implict double integrals over momenta and double sums over frequency; 
the decomposition of the vertex function in terms of a fully two-particle irreducible vertex 
V and T-matrices in particle-hole, crossed particle-hole, and particle-particle channels, 

T = V + Tph + fph + Tpp- (11) 

and a set of Bethe-Salpeter equations relating a T-matrix in the i-th channel to a vertex 
function /j that is two-particle irreducible in that channel 

Ti = kGGh + kGGT,. (12) 
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This set of equations is closed by expressions relating each vertex /j to the T-matrices 
and V. The resulting theory is self-consistent at the level of the single- and two-particle 
propagators. 

We have solved these equations on the CM-5, using a parallel algorithm we describe 
elsewhere (Hess and Serene, 1995), for small Hubbard chains. We employ two addi- 
tional approximations in our parquet approximation calculations that are outside the 
framework of the exact theory: the fully two-particle irreducible vertex V is taken to 
be the bare Hubbard interaction, and a high frequency cutoff is imposed with periodic 
boundary conditions in discrete frequency space (see Deisz, et ai, 1994). Both of these 
approximations may be relaxed (Hess and Serene, 1995). 

In Fig. 6, we show the single-particle propagator G(k, r) calculated in the par- 
quet approximation described above, together with G(k, r) from various single-particle 
self-consistent approximations including the FEA, the shielded potential approximation 
(includes only the particle-hole bubble diagrams together with the second order dia- 
gram), the T-approximation (includes only particle-particle ladder diagrams together 
with the second order diagram), and an approximation including only the second-order 
diagram. The self-consistent calculations were all performed using the same high fre- 
quency cutoff scheme. For comparison we show G(k, r) obtained from quantum Monte 
Carlo (Assaad, 1995). Because of the increased importance of fluctuations in ID, we ex- 
pect that the FEA, shielded potential approximation and T-approximations to be poor. 
Of all the diagramatic approaches shown here, the parquet approximation is closest to 
the quantum Monte Carlo result for a U equal to one quarter of the bare bandwidth. 
In agreement with our sum-rule arguement in §3, the FEA provides a better approxi- 
mation than shielded potential approximation. The approximation containing only the 
second-order 'skeleton' diagram is in better agreement with QMC than any of the other 
single-particle self-consistent theories, at least for this relatively small U. This may in 
part reflect the better agreement with the second moment sum-rule on the single-particle 
spectral function. 

§^ Conclusions 

We have shown that the coupling of spin-fluctuations to quasiparticles, as both evolve 
with decreasing temperature, leads to the formation of a non-Fermi liquid state in the 
2D Hubbard model at half-filling in the fluctuation exchange approximation. The non- 
Fermi liquid state shows a positive slope in the real part of sigma and a local minimum in 
the imaginary part at zero energy. The energy scale over which these anomalies extend 
is roughly the size of the zero temperature gap in mean field theory. We have shown 
that this behavior is preceded by the evolution of a sharp peak the spin-fluctuation T- 
matrix at Q = (tt, tt) and that the evolution of this structure continues with decreasing 
temperature. An analysis of the T-matrix suggests the existence of two temperature 
regimes separated by a characteristic temperature T*. For temperatures above T*, the 
self-energy does not show the anomaly; for temperatures below T*, the anomaly in the 
self-energy is present and the sharp structure in the T-matrix can be described by a 
single temperature dependent energy scale cu* which decreases rapidly with decreasing 
temperature. While the FEA may provide a useful model for coupled quasiparticles and 
spin-fluctuations that self-consistently includes the effects of evolving low-energy anoma- 
lies, it is both desirable and possible to go beyond the FEA. We presented calculations 
that are self-consistent at the levels of both the single- and two-particle propagators for 
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the Hubbard model in ID. Two approximations were employed: replacing the fully two- 
particle irreducible interaction by the bare interaction, and imposing a high frequency 
cutoff. Our results for the single particle propagator suggest the importance of vertex 
corrections in calculating the self-energy and compare favorably with quantum Monte 
Carlo calculations. In future work we will explore the sensitivity of our results to the 
structure of the fully two-particle irreducible interaction and discuss the structure of the 
vertex function. 
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Figure 4: The evolution of the sharply peaked structure of the 
spin fluctuation T-matrix Tgf with decreasing temperature: (top 
panel) The maximum of Tgf at Q and a; = as a function of 
T; (bottom panel) the temperature dependence of the HWHM of 
the sharp peak in Tsf in momentum AqT^/ (❖), the HWHM of 
Tsf{Q,uj) AeTsf (•), and 1/T,/(Q,0) (A). The inset shows the 
merging of ^eTsf with l/Tsf{Q,i 0) below the temperature T* (see 
text) . 
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Figure 5: The fully self-consistent spin response for q = Q = {11,71). Note 
that Xafm{T) does not display the rapid increase with decreasing T that is 
observed in T's/(Q, 0). 
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Figure 6: The single-particle Green's function G(k — 7r/2,r) for an 
eight-site half-filled Hubbard chain ior U — t and T = 0.20 calculated 
in the parquet approximation (solid) compared with the bare propaga- 
tor (dotted), quantum Monte Carlo (A), and the single-particle self- 
consistent perturbation theories: the second order diagram alone (•), 
the T approximation (o), the shielded potential approximation (dot- 
dash), and the fluctuation exchange approximation (dashed). Note that 
of the self-consistent perturbation theory calculations, the parquet ap- 
proximation agrees most closely with the quantum Monte Carlo result. 



14 



